A Conjecture Connected with Units of 
Quadratic Fields 



Nihal Bircan 
December 3, 2012 



Technische Universitat Berlin, Fakultat II, Institut fur Mathematics, MA 
8 — 1 Strasse des 17. Juni 136 D-10623 Berlin, Germany 
bircan@math.tu-berlin.de 

Abstract 

In this article, we consider the order Of = {x + yfy/d : x, j £ Z) 
with conductor / G N in a real quadratic field K — Q(\/d) where d > is 
square-free and d = 2, 3 (mod 4). We obtain numerical information about 
n(f) = n(p) — min{v £ N : e" £ O p } where e > 1 is the fundamental 
unit of K and p is an odd prime. Our numerical results suggest that the 
frequencies of 2 ntp) 01 nfpY should have a limit as the ranges of d and p 
go to infinity. 
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1 Introduction 

Throughout the paper, we consider the real quadratic field K = Q(^/d) where 
d > is square- free and d = 2, 3 (mod 4). An order O in the field if is a subset 
O C K such that 

i) O is a subring of K containing 1 . 

ii) O is a finitely generated Z-module. 
hi) O contains a Q-basis of K. 

Maximal order of an algebraic number field K is the integral closure of the ring 
Z of rational numbers in the field K. So, among various orders of the field K 
there is one maximal one which contains all the other orders which we call Ok- 
Any order O of the field Q(\/rf) is of the form Of = {1, fuj} where / is the index 
[Ok ■ O] namely, the conductor of the order. One can represent any number of 
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the order Of in the form Of = {x + yfVd : x, y € Z}. For d = 2, 3 (mod 4) a 
basis {1, u)} = {1, Vd}. 

For d = 2, 3 (mod 4) one can write very well known Pell's equation as 

(x + yVd)(x — yVd) = ±1 

The s— th solution a; s , y s can be expressed in terms of the first one x±, y\ by 

x s +y s Vd= (xi +yiVd) s (aeN) (1.1) 

The first solution (x%, y{) G N 2 is called fundamental solution to the Pell equa- 
tion. Fundamental units e for real quadratic fields Q(\/d) can be computed 
from the fundamental solution of Pell's equation. So, fundamental units for the 
maximal order Ok of the field K = Q(\/d) are also called fundamental units for 
the algebraic number field K = Q(vd). For more information on units, orders 
and Pell's equation see |Cox[ p. 133 ], |JW1 p. 80 ], |Neul p. 72 ] and other many 
algebraic number theory books. For the efficient computation of units and so- 
lution of Pell's equation see |Len| . One can also refer to |JLW| to see about the 
question how large the fundamental unit e can be. 

In this article, e is the fundamental unit of the real quadratic field Q(y/d) and 
N(e) = ±1 is its norm. We consider the case that the conductor / is an odd 
prime p with p\ d. The aim of this article is to get numerical information about 

n(f) = n(p) = min{v e N : e" e O p }. 

In |BP| . we compute very good upper bounds for n(f) for the case N(e) = +1 
and in |Bir| for the case in particular that N(e) — — 1. It is known that (see 
for instance [BP]) for (^j = =pl, n(f) = n(p) is always a divisor of 2±1 if 

N(e) = +1 and of p ± 1 if N(e) = — 1. Our results become easier to state if we 
consider the integers, quotients q defined by 

q= {£® = +1 where U\ = T l 

\ ffy if*00 = -l where (j) = T l. 

We compute the frequencies of q. Our numerical results suggest that the fre- 
quencies may have a limit as the ranges of p and d go to infinity. We prove a 
theorem suggested by the numerical results and study the expectation of the 
frequency. 

2 Computations and Conjecture 

We recall that we always assume that d = 2, 3 (mod 4) is positive, square-free 
and p is an odd prime. We have Of = {x + yfVd : x, i/gZ}. Every unit of O 
that belongs to Of is also a unit of Of . 
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We compute n(f) = n(p) = min{y e N : e v E Of} for the case that / is an 
odd prime p. For N(e) = +1, the possible maximal values of n(f) are 2±— for 

^ = -1, 2=1 fo r ^ = + i ; These give q = 1 defined in (TQ]| . 

The computer we use has 2x DualCore-Opteron 2218 (2.6 GHz) with 16 GB(1 
Node with 32 GB) RAM. One can also refer to the appendix to see the algorithm. 

Our results depend on two arguments, the range and the quotient q. It would 
be difficult to give a plot. For reasons of space we only present the first 20 
frequencies in our tables as percentages. We take ranges d = [3, . . . , d max \ and 
p = [3, . . . ,p m am] for N(e) = +1. For N(e) = -1, we take d = [2, . . . , d max ] and 
p = [3, . . . ,Pmax\- For N(e) = —1, the computations run faster because there 
are few occurencies so this enables us to implement larger and different ranges. 
To save space and to be neat, we only write d max = Pmax and some of the ranges 
we calculate in the following tables. The Legendre symbol differs according to 
the elements q defined in (|1.2|) . For to e N we define four sequences according 
to the formula (|1.2|) ; 

Fi(g; to) = card{g = (p - l)/2n(p) : d, p < to, (jf) = +1, 7V(e) = +1}. 
F 2 ( g; to) = card{ g = (p + l)/2n(p) : d, p < m, (fj = -1, JV(e) = +1}. 
F 3 ( g ;TO)=card{ ( 7=(p-l)/n(p):d, p < to, (|) = +1, JV(e) = -1}. 

^4(9; m) = card{g = (p + l)/n(p) : d, p < to, (j) = -1, iV(e) = -1}. 

For instance, for the case iV(e) = +1, by the definition (|1.2[) of g, we have 
q < 2^—. Hence Fi = for <j > . This means that, when the ranges become 
larger, more and more frequencies will become positive. This suggests that for 
a fixed q the frequencies will slowly fluctuate. 

Furthermore, let Sj(m) — Fj(q]m). In the tables values shows the number 
of occurencies of (^j = +1 and (^j = — 1 respectively. The tables suggest the 
following conjecture; 

Conjecture 2.1. Let j = 1,2,3,4. There is a probability distribution Pj(q) 
such that, for all q, 

Fj{q; m)/Sj(m) — > Pj(q) as to — > 00. 
More precisely: For j = 1, . . . , 4 there is a function 

Pj ■ : N -> [0, 1] w'ift P A<1) = 1 

such that for every 5 > £/iere exists Too wii/i tte property that 

\F j (q;m)/S j (m)-P j (q)\<S 

for all q and m > mo . 
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iV( £ )=+l,(|) 


= +1 


i 


10000 


20000 


40000 


60000 


80000 


100000 


1 


57.3 


56.9 


56.9 


56.7 


56.7 


56.6 


2 


11.8 


12.0 


12.0 


12.1 


12.2 


12.2 


3 


9.91 


9.89 


9.89 


9.88 


9.86 


9.88 


4 


4.66 


4.72 


4.64 


4.64 


4.65 


4.67 


5 


2.86 


2.84 


2.84 


2.84 


2.82 


2.82 


6 


1.98 


2.05 


2.10 


2.13 


2.14 


2.16 


7 


1.34 


1.35 


1.34 


1.35 


1.34 


1.34 


8 


1.12 


1.16 


1.16 


1.16 


1.16 


1.16 


9 


1.09 


1.08 


1.09 


1.08 


1.09 


1.09 


10 


0.604 


0.610 


0.609 


0.612 


0.618 


0.619 


11 


0.530 


0.515 


0.521 


0.518 


0.517 


0.510 


12 


0.817 


0.823 


0.841 


0.854 


0.866 


0.874 


13 


0.349 


0.365 


0.363 


0.362 


0.363 


0.363 


14 


0.294 


0.292 


0.289 


0.292 


0.291 


0.295 


15 


0.497 


0.489 


0.494 


0.492 


0.491 


0.495 


16 


0.280 


0.306 


0.302 


0.296 


0.293 


0.293 


17 


0.216 


0.214 


0.210 


0.214 


0.213 


0.210 


18 


0.221 


0.229 


0.237 


0.239 


0.238 


0.240 


19 


0.161 


0.164 


0.159 


0.164 


0.163 


0.162 


20 


0.233 


0.238 


0.239 


0.236 


0.238 


0.240 


Values 


2249621 


8330759 


31140582 


67496484 


116576513 


178609439 



N(e) = +!,(!) 






10000 


20000 


40000 


60000 


80000 


100000 


l 


56.3 


56.2 


56.0 


56.1 


56.1 


56.1 


2 


14.5 


14.4 


14.4 


14.3 


14.3 


14.4 


3 


9.94 


9.97 


9.98 


9.97 


10.0 


9.98 


4 


3.32 


3.31 


3.36 


3.37 


3.37 


3.34 


5 


2.77 


2.76 


2.80 


2.80 


2.80 


2.81 


6 


2.45 


2.48 


2.49 


2.49 


2.49 


2.50 


7 


1.32 


1.36 


1.35 


1.36 


1.36 


1.35 


8 


0.868 


0.840 


0.849 


0.836 


0.839 


0.837 


9 


1.07 


1.10 


1.10 


1.11 


1.12 


1.12 


10 


0.718 


0.728 


0.721 


0.725 


0.736 


0.729 


11 


0.512 


0.502 


0.518 


0.518 


0.519 


0.514 


12 


0.561 


0.566 


0.583 


0.581 


0.578 


0.571 


13 


0.352 


0.363 


0.352 


0.355 


0.351 


0.354 


14 


0.341 


0.339 


0.328 


0.339 


0.333 


0.332 


15 


0.480 


0.487 


0.508 


0.507 


0.508 


0.506 


16 


0.202 


0.200 


0.206 


0.206 


0.210 


0.209 


17 


0.190 


0.199 


0.201 


0.201 


0.200 


0.204 


18 


0.274 


0.268 


0.274 


0.274 


0.272 


0.273 


19 


0.165 


0.168 


0.169 


0.168 


0.168 


0.165 


20 


0.177 


0.176 


0.175 


0.175 


0.174 


0.171 


Values 


2272057 


8390244 


31303879 


67789746 


117013651 


179198341 
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N(e) 


= - 1 »(l) 


= +1 






i ^^^^ 


40000 


60000 


80000 


100000 


200000 


400000 


1 


37.8 


37.6 


37.6 


37.5 


37.4 


37.5 


2 


18.7 


18.8 


18.8 


18.8 


18.7 


18.7 


3 


6.58 


6.56 


6.60 


6.62 


6.64 


6.62 


4 


14.0 


14.0 


14.0 


14.0 


14.0 


14.0 


5 


1.90 


1.91 


1.89 


1.90 


1.89 


1.90 


6 


3.28 


3.30 


3.29 


3.29 


3.31 


3.31 


7 


0.891 


0.897 


0.899 


0.898 


0.894 


0.892 


8 


3.46 


3.48 


3.49 


3.51 


3.50 


3.49 


9 


0.724 


0.716 


0.725 


0.728 


0.733 


0.735 


10 


0.930 


0.928 


0.936 


0.937 


0.943 


0.938 


11 


0.352 


0.351 


0.348 


0.340 


0.343 


0.339 


12 


2.49 


2.50 


2.48 


2.48 


2.49 


2.49 


13 


0.237 


0.239 


0.240 


0.243 


0.239 


0.245 


14 


0.440 


0.440 


0.438 


0.442 


0.444 


0.447 


15 


0.329 


0.329 


0.330 


0.333 


0.333 


0.334 


16 


0.866 


0.868 


0.863 


0.866 


0.872 


0.874 


17 


0.140 


0.143 


0.143 


0.140 


0.138 


0.140 


18 


0.371 


0.371 


0.365 


0.366 


0.369 


0.368 


19 


0.105 


0.109 


0.108 


0.107 


0.109 


0.109 


20 


0.705 


0.706 


0.705 


0.702 


0.706 


0.703 


Values 


2812857 


5965306 


10163034 


15404441 


56042335 


205444859 







N(e) 


— • (i) 


= -1 








40000 


60000 


80000 


100000 


200000 


400000 


l 


37.9 


37.6 


37.6 


37.6 


37.6 


37.5 


2 


37.1 


37.2 


37.2 


37.2 


37.3 


37.3 


3 


6.66 


6.65 


6.64 


6.62 


6.61 


6.65 


4 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


5 


1.90 


1.90 


1.92 


1.90 


1.90 


1.89 


6 


6.59 


6.60 


6.64 


6.63 


6.64 


6.64 


7 


0.889 


0.910 


0.897 


0.891 


0.891 


0.899 


8 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


9 


0.727 


0.732 


0.729 


0.723 


0.731 


0.733 


10 


1.85 


1.85 


1.84 


1.87 


1.87 


1.88 


11 


0.344 


0.345 


0.347 


0.351 


0.338 


0.343 


12 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


13 


0.233 


0.236 


0.238 


0.240 


0.241 


0.239 


14 


0.895 


0.907 


0.907 


0.900 


0.899 


0.895 


15 


0.342 


0.336 


0.336 


0.332 


0.334 


0.333 


16 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


17 


0.135 


0.138 


0.139 


0.140 


0.138 


0.138 


18 


0.734 


0.741 


0.750 


0.754 


0.742 


0.742 


19 


0.112 


0.112 


0.113 


0.111 


0.110 


0.110 


20 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


Values 


2828439 


5992331 


10206629 


15464072 


56197833 


205855014 



Comments: The first two tables appear to be rather similar. In the first two 
tables q = 1 predominates. In the third table the sum of frequencies for q = 1 
and q = 2 is about the frequency for q = 1 of the first two tables. 

But the last table is quite different. The first two frequencies are almost equal 
and we notice that when q = 4fc for k € N, ^(4; m) — 0. Now. we give a proof 
that this is indeed true. 



Theorem 2.1. Let N(e) = — 1, (J J = — 1. T/ien £/ie case q = 4 does noi 
appear. 

Proof. Suppose that a = 4 does occur. Then i*4(4;m) ^ and thus v = 
In |Bir| we wrote 

e" = A„ + B^vVd 
In our case v — 2±- and thus Bp+i _ r Since !/ e N it follows that p = 3 
(mod 4). In |Bir| Theorem 3.1 we proved that, for (^j = —1, Ap±i = (mod p) 
if (jjfj = — 1 ail d for iV € Z, TV 7^ by using the following formula 

2 (x 2 - N) B v ^{xf = A 2v {x) - N» 

we obtain 

2 (x 2 + l) Bj+i^fi) 2 = As+iix) - (— l) 2 ^ 
= 1 (mod p) 

So, ^ (mod p) which is a contradiction. □ 

4 

We compute the expectation of the probability distribution for N(e) = +1, 
= +1 which are listed below. 



Ranges= m 


1000 


10000 


15000 


20000 


25000 


30000 


Expectation= E 


3.921 


6.086 


6.397 


6.759 


6.946 


7.098 


log(-m) X log(logm) 


13.344 


20.448 


21.758 


22.697 


23.441 


24.047 


E/(\og(m) X log(logm)) 


0.293 


0.297 


0.294 


0.297 


0.296 


0.295 



Ranges= m 


40000 


50000 


60000 


70000 


80000 


100000 


Expectation= E 


7.300 


7.472 


7.681 


7.785 


7.922 


8.152 


log(m) X log(logm) 


25.012 


25.763 


26.383 


26.536 


27.362 


28.129 


E/(log(m) X log(logm)) 


0.291 


0.290 


0.291 


0.293 


0.289 


0.289 



This table suggests that the expectation has the order of log(7«) x log(logm) 
for large m where m is the range of d and p. 
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Appendix 



We use the following procedure to obtain the values. The computation is similar 
for N(e) = — 1. At each step of the computation the numbers are modulo p so 
that they are 32 bit numbers. To speed things up, one should not calculate e l 
as an element of the number field. It is faster to just calculate the coefficients 
modulo p with respect to the basis {1, Vd}. They can be calculated with a fast 
exponentiation algorithm. Computations are made by using |BCP] . 

Algorithm 1 Finding n(f) = n(p) 



1: Compute the Norm of e and the Legendre Symbol Is 4— l^j 

2: if Norm(e) = 1 then 

3: q <- (p - ls)/2. 

4: else 

5: q <— p— Is 

6: end if 

7: Calculate factorization of q — p^ . . .p\* and set v <— 1 
8: for i in [l..f] do 

9: Calculate F[i) max{6 £ {0, r^}| for which the second coefficient 
y q / p b (see (| 1 . lj) ) of e pb i is divisible by p} 

F[i] 
V v ■ p i 

10: end for 

11: Compute for N(s) = lo r £| for N(s) = -1. 
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